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Abstract 

Riemannian geometrical tools, such as Ricci collineations and Killing symmetries , so often 
used in Einsteins general theory of gravitation are here applied to plasma physics to build 
magnetic surfaces from Einstein plasma metrics used in tokamak devices. It is shown that the 
Killing symmetries are constrains the Einstein magnetic surfaces while the Killing vectors are 
built in terms of the displacement of the toroidal surface. The pressure is computed by applying 
these constraints to the pressure equations in tokamaks. A method, based on the sign of the 
only nontrivial constant Riemann curvature component, is suggested to classify tearing insta- 
bility. Throughout the computations two approximations are considered: The first is the small 
toroidality and the other is the small displacement of the magnetic surfaces as Einstein spaces. 
PACS numbers: 02.40.Hw-Riemannian geometries 
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I Introduction 



The tools of the Riemann geometry [1] so often used in other important areas of physics, such 
as Einstein theory of gravitation, have also been used by Mikhailovskii [2] to investigate the 
tearing and other sort of instabilities in confined plasmas [3], where the Riemann metric tensor 
played a dynamical role interacting with the magnetic field through the magnetohydrodynam- 
ical equations (MHD). More recentely Garcia de Andrade [4] has applied Riemann metric to 
investigate magnetic flux tubes in superconducting plasmas. Earlier Thiffault and Boozer [5] 
have investigated the Riemann geometry in the context of chaotic flows and fast dynamos. In 
this paper we use the tools of Riemannian geometry, also used in other branches of physics 
as general relativity [6], such as Killing symmetries , Riemann and Ricci [7] collineations , 
shall be applied here to generate magnetic nested surfaces in Tokamaks through the built of 
Einstein spaces obtained from Tokamak plasma metric. This work is motivated by the fact 
that the magnetic surfaces are severely constrained in tokamaks [3] and Killing symmetries are 
well applied every time we have symmetries in the problem as in solutions of Einstein equa- 
tions of general relativity. Equilibrium of these surfaces or their instabilities are fundamental 
in the constructio of tokamaks and other plasma devices such as stellarators where torsion is 
also present. The magnetic surfaces are more easily obtained when symmetries are present. 
This is our main motivation to apply the special Riemann geometrical techniques of Killing 
symmetries and Ricci collineation to obtain magnetic surfaces formed by Einstein spaces. To 
simplify matters we shall consider two usual approximations from plasma physics [3] which are 
the small toroidality or inverse aspect ratio e = -| << 1 where here R represents the external 
radius of the torus and a is its internal radius, and the Shafranov displacement A << 1 as 
well. The pressure of the tokamak is obtained from the tokamak Shafranov shift equation. 
Constant pressure closed ergodic nested surfaces in magnetohydrostatics have also been shown 
by Schief [8] to be generated by solitons. This is another mathematical technique , distinct 
from ours, is another use of mathematical theory to generate of nested magnetic surfaces in 
plasmas. The paper is organised as follows: In section 2 we review the Riemannian tech- 
nique of Killing vector and Ricci and Riemann colineations not usually familiar to the plasma 
physicists and solve the Ricci tensor components from the plasma metric by considering that 
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nested surfaces are formed by Einstein spaces, where the Ricci tensor is proportional to the 
metric, and solve the Ricci collineation equations to find out the Killing vectors stablishing a 
geometrical method for the classification of tearing instabilities. Conclusions are presented in 
section 3. 



II Ricci collineations from plasma metrics and tear- 
ing instabilities 

dynamos. Let us now start by considering the plasma metric given by Zakharov and Shafra- 
nov [9] to investigate the evolution of equilibrium of toroidal plasmas. The components g ik 
(i,k=l,2,3) and (a,9,z) as coordinates, of their plasma metric are 

flu = 1 - 2A'cos6 + A' 2 (II.l) 

922 = a 2 (II.2) 

^33 = (R - A + acosO f (II.3) 

0i 2 = aA'sinO (II.4) 

where z = asinO and the dash represents derivation with respect to a. In our approximation 
the last term in the expression (II.l) may be dropped. The magnetic surface equations are 
tori with circular cross-section and equations 

r = R- A(a) + acosO (II.5) 

z = asinO (H-6) 

Let us now compute the Riemann space of constant curvature represented by the Riemann 
tensor components 

Ri jk i = A(g ik gji - gug jk ) (H.7) 

where A ia constant which is called de Sitter cosmological constant. Contraction of expression 
(11.7) in two non-consecutive indices, otherwise the symmetry of the Riemann curvature tensor 
Rijki = —RjiH = Rjiik would make them vanish, yields the Einstein space Ricci relation 

Rik = 2Ag tk (11.8) 



The Ricci collineations equations are given by 

[m^rf + R a d k rf + R kl d iV l = (II.9) 

where rf are the components of the Killing vector fj which defines the symmetries of the 
associated space, and di := are the components of the partial derivative operator. This 
equation is obtained from the more elegant definition in terms of the Lie derivative C v as 

C v Rik = (11.10) 

In the next section we shall solve the Ricci collineation equations in terms of the plasma metric 
above. 



Ill Nested surfaces as Einstein spaces in plasmas 

Let us now consider the application of the above plasma metric into the Ricci collineation 
equation, which yields 

[d l R 11 ]rf + 2R u d 1 rf = Q (III.ll) 

[d l R 22 } V l + 2R 2l d 2 r ] l = (111.12) 

We shall consider now just two independent coordinates (x 1 = a, x 2 = 0) since nested surfaces 
are bidimensional in the case of plasmas, Let us now compute the Riemann tensor components 
in the linear approximation 

Rl212 = | ^912 _ d 2 g u _ d 2 g 22 
dx 1 dx 2 dx 2 dx 2 dx 1 dx 1 

Substitution of the plasma metric above into expression (111.13) yields the expression 

R 1212 = [2 + (asinO + 2cos6)} (III. 14) 

It is easy to show that the components -R1313 and i? 2 323 both vanishes within our approx- 
imations. At this point we consider that 8 is so small that sin8 vanishes and cos8 = 1 
this simplifies extremely our metric and turns it into a diagonal metric where g± 2 = and 
g bb — (gtb)^ 1 (b = 1,2) and this allows us to compute the components of the Ricci tensor 
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from the Riemann component. But before that let us compute the use the condition that the 
nested surface is an Einstein space to compute the Riemann component again 

R 1212 = Aa 2 [1 - 2 A'cosO] (III. 15) 

Since both expressions for the Riemann component -R1212 must coincide, equating expressions 
(111.15) and (111.14) yields an expression for the derivative of the Shafranov shift A as 

A' = [l + ^] (111.16) 
Integration of this expression yields the value of the shift in terms of the radius a as 

A = [a+^] (III. 17) 

Since the radius a is assumed to be small we may neglect terms of the order a 3 which yields 

A = a (111.18) 

which satisfies the well-known boundary condition A(0) = 0. From these expressions one may 
also compute A" = Aa. An important result in plasma physics is that tearing instabilities 
coming from ion or electron currents possess the shift condition A' < 0. This condition would 
be clearly violated from expression (111.16) unless the A curvature constant would be negative 
and in modulus > 1. Note that this situation is very similar to the condition of favorable 
or unfavorable curvature for the instabilities in plasmas [3]. The main difference is that here 
we are referreing to Riemann curvature and not to Frenet curvature of the magnetic lines in 
plasmas. This suggests another method to classify geometrically tearing instabilities. Actually, 
since has been shown [11] recently that the Riemann tensor in plasmas can be expressed in 
terms of the Frenet curvature both methods seems to be equivalent. Now let us compute the 
Ricci components Rn and .R22 from the component -R1212 by tensor contraction with metric 
components g 11 and g 22 , which results in the expressions 

Rn = — + A'] (111.19) 

cr 

and 

R 22 = [2 + 6A'] (111.20) 



which in turn yields the expressions 



diRn = ~[1 + A'] - A[A»] = A (111.21) 



a 3 a 2 a 3 

and 

fibril = -— = - l 2 J (111.22) 
a a 

^1^22 = A' (111.23) 

d 2 R 22 = aA' (111.24) 

Substitution of these derivatives of the Ricci tensor components into the Ricci collineations 
equations one obtains the following set of PDE equations 

[a- 1 + j}d lV 1 -r] 1 ^ = (111.25) 

which yields 

rf = exp[-} (111.26) 
a 

If} i 

—d 2 rf + rf + 6Aaexp\-} = (111.27) 
a a 

which by considering the gauge d 2 rf = one obtains the constraint 

rf = -QA exp [-] (111.28) 

CL 

for the Killing vector field poloidal component. The relation 

\ = --A' 1 (111.29) 
rf 6 V ; 

allows us finally to write down the Killing vector as 

ff=-]:A- 1 [l,-6A]Ti e (111.30) 
o 

which depends again upon the curvature scalar A which is central in the classification of tearing 
instabilities. 
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IV Conclusions 



In conclusion, we have investigated a method of classification and identification of tearing 
instability, allowing for example to distinguish between tearing instabilities that comes from 
ions and electron currents or not, based on the Riemann curvature constant submanifolds as 
nested surfaces in Einstein spaces. The Killing symmetries are shown also to be very useful 
in the classification of plasma metrics in the same way they were useful in classifying general 
relativistic solutions of Einstein's gravitational equations in four-dimensional spacetime [6]. 
Other interesting examples of the utility of th is method is the Ricci collineations investigations 
of the twisted magnetic flux tubes and the Arnolds metric for the fast dynamo [10, 11, 12]. 
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